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Abstract 



^ ■ We propose a new lattice action for non-abelian gauge theories, 

, which will reduce short-range lattice artifacts in the computation of 

■ the topological susceptibility. The standard Wilson action is replaced 

by the Wilson action of a gauge covariant interpolation of the original 

■ fields to a finer lattice. If the latter is fine enough, the action of all 

configurations with non-zero topological charge will satisfy the contin- 
uum bound. As a simpler example we consider the 0(3) a-model in 
two dimensions, where a numerical analysis of discretized continuum 
instantons indicates that a finer lattice with half the lattice spacing of 
the original is enough to satisfy the continuum bound. 
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Field configurations with non-zero topological charge are expected to 
have a strong influence on the dynamics of asymptotically free theories. In 
QCD, such configurations are responsible for breaking axial symmetry and 
resolving the U(l) problem The study of these effects however requires 
non-perturbative techniques and one would expect that ultimately Monte 
Carlo methods on the lattice would be best suited to it. The observable 
to consider is inspired by the classic large-iV c analyses of Witten and of 
Veneziano Q, which showed that, 

2 i 2 r, 2 6 Xt / -i \ 

ny + m v - 2 m K = — T , (1) 

Jit 

where xt is the topological susceptibility. In the continuum it is given by, 



Xt= d X < q(x)q(0) > \no quarks (2) 



with q(x) being the topological charge density, 



<?0) = — ^ W^ Tr I F ^ F p«]- ( 3 ) 



The topological charge, Q = J q(x), is an integer if the field strength vanishes 
at infinity or if (euclidean) space-time is compact. A continuum analysis also 
shows that the action of any configuration with non-zero topological charge 
must satisfy the following bound, 

% 

A big effort has been devoted to the study of the topological susceptibil- 
ity on the lattice. There are several choices for the operator q(x). The naive 
discretization of (]3|) does not yield integer values for Q and requires renor- 
malization factors ||. On the other hand, the cleaner geometrical definition 
due to Liischer [Q] gives an integer-valued topological charge and does not 
require renormalization. Here we will deal only with a geometrical definition 
very similar to Liischer's. The geometrical topological susceptibility is then 
obtained by computing, 

<Q 2 > ~ 

Xt = y , 

where V is the volume of the lattice. This definition is clearly equivalent to 
(Q) in finite volume. From the continuum formula (|l]), it follows that the 
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topological susceptibility in QCD should scale as (mass) 4 in the continuum 
limit, namely 



Xt 



' l exp ) 4 , (6) 



where b is the lattice spacing and (3 (go) ^ — AjSo is the leading term of the 
beta function. 

However, it was found, first in the 0(3) cr-model in two space-time di- 
mensions H pi and then in non-abelian gauge theories in four jjj] Q , that the 
standard actions in both cases give rise to short-range fluctuations with non- 
zero geometrical topological charge and smaller action than the continuum 
bound. These fluctuations, often referred to as "dislocations", overwhelm 
the contribution of the slowly varying fields, which would otherwise dom- 
inate in the continuum limit, and are expected to destroy the scaling of 
the topological susceptibility. On the other hand, if the bound is satisfied, 
the semiclassical continuum analysis of non-abelian gauge theories indicates 
that the susceptibility should indeed be ultraviolet finite, and therefore scale. 
Satisfying the continuum bound is thus a sufficient condition for scaling in 
non-abelian gauge theories ||] @ ||] f\ This suggests that eq. (||) can be 
satisfied on the lattice by giving dislocations a larger action. 

Several proposals to solve this problem have been considered in the lit- 
erature, like the cooling method || and the use of improved actions jjj [1C]. 



In this paper we propose to use a new action which satisfies the continuum 
bound. It is related in spirit to the actions proposed in [TfJ, but may be 
simpler to implement. 

The idea behind our new action is easy to understand once one realizes 
the mismatch between the geometrical definition of topological charge 
and the Wilson action for gauge fields which allows dislocations to arise. 
Consider a continuum instanton A^(y) which saturates the bound (Q), and 
discretize it on a lattice of spacing b, 

rsb+bft 

Up(s) = Pexp(i / dyA^y)), (7) 

J sb 

where s b are the sites of the b lattice. The geometrical definition of topo- 
logical charge assigns a non-zero value even to a lattice configuration (^) 

*It has been argued by Liischer|| that in the case of the 0(3) a-model, the failure in 
finding scaling is not related to the existence of dislocations and is an essential problem of 
this model: the partition functional in the instanton sector shows a logarithmic ultraviolet 
divergence in the continuum analysis. 
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obtained from very small instantons, of order of the lattice spacing b. On 
the other hand, it is clear that the Wilson action very poorly approximates 
the continuum action of the original instanton of 0(b), and in fact it is 
generically smaller and therefore does not satisfy the bound. Such rough 
configurations are dislocations and can destroy the scaling of the suscepti- 
bility. An important observation is that the geometrical topological charge 
assigned to a lattice configuration can be understood as the naive topological 
charge of a continuum configuration obtained by smoothly interpolating the 
lattice configuration. Then it is clear that if, instead of using the standard 
Wilson action of the original lattice configuration, we use the continuum 
action of the interpolated configuration, the continuum bound is necessarily 
satisfied, as first suggested in [|Tl|], 

More concretely, in [^] we described a procedure to obtain a continuum 
gauge field a fJi (y) which interpolates any 6-lattice configuration^]. The inter- 
polation is local and gauge covariant, i.e. for a 6-lattice gauge transformation 
fi(s), 

a,[U n }=a^[U], (8) 

where w is a gauge transformation in the continuum. (Other interpolations 
can also be used in this context [PH-) A geometrical topological charge of 
the 6-lattice configuration is defined as the one associated to the interpolated 



field 112 



This definition has the same properties as Liischer's original definition B 
and requires roughly the same computational effort. Now, it is clear that 
replacing the standard Wilson action by the continuum action of a^, i.e. 

s cont = -L( d % Tr[f^(y)f^(y)], (10) 

insures that the continuum bound is satisfied for the same reason as it is in 
the continuum. Notice that this is a perfectly gauge invariant action for the 
lattice field t/ M (s), by eq. @. 

^The continuum field is differential) le inside each 6-hypercube and transversely con- 
tinous across the 6-boundaries, and when discretized according to (Q) gives back U l _ l {s) 

& 
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However, for numerical purposes the continuum action is impractical 
LlJ. The central observation of this paper is that it is not necessary to 
interpolate all the way to the continuum, but just to a finer lattice, with 
lattice spacing / (we will take bj f to be integer). From now on, we refer to x 
as the sites on the /-lattice and s as the sites on the ^-lattice ( x^ = s^+m^^, 
= 0, 4 — 1). The interpolation procedure in |l2| gives a set of link 
variables u fJl [U](x), such that 

- - log(u^x)) = a^x) + 0(//6 2 ), (11) 

where a^{x) is the continuum interpolation discussed above, at point x. 
On the / lattice, we can simply choose the standard Wilson action. The 
partition functional will then have the form, 

Z = fl[VU(s) e-^I*] (12) 

•* s 

where, DU is the usual Haar measure for non-abelian gauge fields on the 
6-lattice, and the Wilson action in terms of the interpolated link variables 
u[U] is given by, 

sLsonHu]] = 4 E E i 1 - *vN + h.c), 

Ufj, u [x] = Tr[u^(x)u u (x + fj)u^(x + u)ul(x)]. (13) 

Again, this action is gauge invariant, because the functional u[U] is gauge 
covariant [12|. From eq. (|TT|) it then follows that, 



sLson = y / A Tr\SMSM\ + 0(f/b). (14) 

Although the O(ffb) terms are not necessarily positive definite, clearly by 
taking f/b small enough we can come arbitrarily close to satisfying the 
continuum bound, so that Xt shows scaling. Determining how small this 
ratio must be in practice requires a numerical analysis, which is beyond the 
scope of this letter. However, we have carried out a numerical analysis of this 
issue in a simplified model in two dimensions, which shares many features 
with four dimensional Yang-Mills theories and our results encourage us to 
believe that t need not be very small in order to recover scaling of \t- 



4 



1 0(3) d-Model in 2D 



The 0(3) er-model in two dimensions [O] |14] is the simplest asymptoti- 
cally free field theory and, as is well-known, it has instanton solutions [13]. 
A continuum bound on the action exists for configurations with non-zero 
topological charge. We will show that, while in the standard lattice for- 
mulation of this model ||, there are topologically non-trivial configurations 
with a smaller action than the continuum bound, this is not the case when 
we use an improved action along the lines described above. For a similar 
discussion in the context of renormalization improved actions see [15]. We 



will not address in this paper the interesting question of whether short- 
range configurations satisfying the continuum bound in this model can still 
dominate in the continuum limit, as argued by Liischer Q. 

The 0(3) o"-model in the continuum is defined by the action, 

S = A fd 2 xJ2(d»N(x)) 2 , (15) 

where N is a 3-component real field satisfying the constraint N =1 (which 
defines the surface of a sphere of unit radius). The continuum topological 
charge in this model is given by, 

Q = ^J d2x ^ ■ x &^)> (16) 

which is the number of times that space-time wraps around the iV-sphere. 
Using the identity [ffq], 



j(8^N + e^N x d v Nf = i(c^iV) 2 - X - e^N • d^N x d u N, (17) 



it follows that, 

S=^ + A f d 2 x(d^N + e^[Nxd u N}) 2 , (18) 
and, since the second term is positive definite, 

S>^. (19) 

% 



Continuum instantons saturate this bound [ 13 1 
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In a standard lattice treatment the action is, 



s b = AEE(W)) 2 (20) 



S fl 



where s are the sites of a two dimensional lattice and d^N = N(s+fi) — N(s). 
The topological charge as defined by Berg and Liischer || is given by, 

Q 6 = E*). (21) 



f[si,s 2 ,s 3 



^-{Im[Log[f[s, s + p,,s + pL + v) f[s, s + fi + v,s + £>]]]}, 
1 + E i<3 N( Si ) • N( Sj ) + iJV(si) • [JV(s 2 ) x i\f(a 3 )] 



[2 n 4<J (i + ^)-^-))] 1/2 

(22) 



It is not hard to understand this formula. Consider the plaquette (s, 1,2). 
The spin variables at the corners are four points on the sphere. We can form 
two triangles with corners at these points and with sides along geodesies, 
Ti(s) = {N(s),N(s+i),N(s + l + 2)) andT 2 (s) = (N(s), N(s + l + 2), N(s + 
2)). q(s) is simply the sum of the area of these two image triangles on the 
iV-sphere. If periodic boundary conditions are imposed, Q b is necessarily an 
integer. 

With these definitions of the action (]20| ) and the topological charge ( |22|) , 
the analysis in J5||| showed that there are dislocations, i.e. configurations 
with non-zero topological charge that have a smaller action than the contin- 



uum bound (19). We will show that this picture changes considerably when 
the action of an interpolation of TV" is used. 

We first interpolate the lattice spin variables to a finer lattice, by moving 
along geodesies on the sphere. The interpolation will be done locally, i.e. 
the interpolated fields within a plaquette (s, 1,2) of the b lattice will only 
depend on the four spins associated with this plaquette, i.e. N(s),N(s + 
l),N(s + 2), N(s + 1 + 2). The points of the finer lattice contained in the 
plaquette at s, are x = s + til + i 2 2, with < t\p < 1 and multiples of f /b. 
The spin fields at these points will be denoted by n(ti,£ 2 )- We first define 
the interpolated spin fields at the corners in the obvious way, 

n(t 1 ,t 2 ) = N(s + t 1 l + t 2 2) ti,t 2 = 0,l. (23) 



6 



We now interpolate along the one dimensional boundaries of the plaquette, 
by moving along geodesies on the sphere. For < t\ < 1, 

m sin[fl(0,0;l,0)(l-t 1 )] ^ n sin[fl(0, 0; 1, 0) h] ^ 
n(tl ' 0) = sin[6>(0,0;l,0)] n( °' 0) + sin[f9(0,0; 1,0)] " (1 ' 0) 



n(t1 ' 1} = 8in[(9(0, 1; 1, 1)] n( °' 1} + sin[fl(0, !;!,!)] w(1 ' ^ 

(24) 

where we have defined, 

9(t 1 ,t 2 ;t' 1 ,t' 2 ) = d [n(ti,*2),n(t / i,t , 2 )] = arccos[n(ti,* 2 ) ■ n(0' 2 )]- (25) 

The expression for spins on the boundaries along direction 2 are analogous 
so we skip the formulae. 

We now define the spins in the interior. In order to do this we first 
interpolate along the diagonal geodesic connecting N{s) and N(s + 1 + 2), 

_ S m[*(0 0;l,l)(l-<)] BinP(0,0;l,l)t] 
V ; sin[0(O,O;l,l)] ^ ' ' sin[0(O, 0; 1, 1)] V ' ; V ' 

For the interior points with t\ > t 2 , we obtain the spin by simply moving 
a fraction t 2 of the distance along the geodesic from n(t\,0) to n(t\,t\). 
Similarly for points with t 2 > ti, the spin variable is obtained by moving a 
distance t 2 along the geodesic linking n(ti, 1) and n(t\, t\). This corresponds 
to separately interpolating the interior points of the two image triangles 
Ti(s) and T 2 (s). Defining, 

a(ti) = 0(ti,O;ii,ti) /S(ti) = e(ti,*i;ti,l), (27) 

the final expression is, 

sin [a (1 — t 2 /U)} _, . sin [a *2/*il 
sm|a] sm|a] 

act..*,) = sm[g " i) + *.>*. P8) 

sin[p] sm[p] 

with ti = 1 — U. 



7 



The interpolation ( p8| ) is now a configuration on an /-lattice and we 
define the improved action to be, 

S J = ^EE («x)) 2 , (29) 

Z % x ^=1,2 

where x is any point on the /-lattice and <9^n(x) = n(x + / /bfi) — n(x) = 
0(f/b). 

Now, we can easily derive a bound on the action, noticing that by con- 
struction, 

Q f = Q\ (30) 

where Q$ is the topological charge for the /-lattice configuration: it is given 
by (|2~2"D, simply substituting N by n and s by x. Now we can expand in 
d n ~ 0(/) and find that, 

Q J = j- E "(a?) • x 4n(x)) + 0(f/b). (31) 

Using ([H]) one finds, 

S f = + A E + e /4™ x <in]) 2 + 0(//6). (32) 

Then, 

S f>^. + 0{f/b). (33) 
So 

Since it is not possible to prove that the 0(f /b) terms are positive definite, 
we do not know how small an f/b we need to be sufficiently close to the 
continuum bound. 

In order to address this question, we considered the discretization of 
a continuum instanton configuration with unit topological charge given by 



m + in 2 = 2w/(l + \w\ 2 ) n 3 = (1 - \w\ 2 )/(l + \w\ 2 ) 



w{xi + ix 2 ) = — — (34) 

Xl + — C(l + 2 J 
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where n = (ni,ri2, W3). We take a and c to be real and define, 

a = (r + r/2\/2) c = (r - r/2\/2) (35) 

The radius of the instanton is proportional to r = \a — c\, while the center is 
located at (rQ,rQ), with ro = \a + c\/2. Generically there is always a critical, 
r c , below which Q b = 0. We expect, and find numerically, that r c ~ b. 
We consider the center to be situated at the center of the volume to reduce 
finite volume effects. Obviously, the continuum instanton is not periodic. 
We impose periodic boundary conditions by defining, 

\ — -r \ 1L max sin(7r(l — Si L max )) 

N( Sl ,s 2 ) = n{z 1 ,z 2 ), Zi = — -, 36 

vr 1 - cos(7r(l - Si/L max )) 

where Z{ are coordinates on the infinite plane and Sj G (—L max ,L max ) are 
the coordinates on the lattice (torus). (2L max ) 2 is then the number of lattice 
sites. The connection between these two sets of coordinates is established 
through sterographic projection. This deformation of the infinite volume 
instanton is small if the instanton size is much smaller than L max . For such 
configurations therefore, the continuum action nearly saturates the bound. 
An alternative procedure would be to discretize the continuum instanton 
solutions of this model in a torus [15]. 



Figure 1 sumarizes our results. It represents the action of the discretized 
instanton configuration as the radius is varied. The continuous line corre- 
sponds to the standard action (p0|), while the dashed lines correspond to 
the improved action for different values of the ratio / jb. It is clear that the 
standard action is problematic, since for r > r c the action is smaller than 
the continuum bound. For the new action however the situation is differ- 
ent. Not only is the continuum bound satisfied, but also as the instanton is 
shrunk to sizes of 0(b), a small barrier develops, separating the Q b = and 
Q b = 1 sectors. The existence of this barrier is easy to understand. The new 
action is the action of the instanton which has been first discretized on the 
6-lattice and then interpolated. For a large instanton (compared to b), the 
interpolation should recover approximately the original configuration and 
so the action must be near that of a continuum instanton and show scale 
invariance. On the other hand, for an instanton of size ~ b a lot of infor- 
mation is lost in the discretization and the interpolation is not expected to 
give a configuration similar to the original continuum instanton. In general, 
the interpolation of small discretized instantons of 0(b) will then be some 
other configuration that need not be even approximately an instanton, and 
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cy qQ =1 

^ I '-'min 



1.2 



. 8 



. 6 



0.4 



0.2 




r/6 



Figure 1: Action of a discretized instanton of the 0(3) model (normalized 
to the continuum one-instanton bound ([19|)) as a function of its radius, r. 
The full line is the standard action in a 100 x 100 lattice and the dashed 
lines correspond respectively to f/b = 1/2,1/4,1/6 (smaller f/b, smaller 
dashing). The vertical line at r c ~ 1.46 separates the Q b = 0, 1 sectors. 



consequently its action will be larger than the bound, since only instantons 
saturate the bound. Thus the action must increase as we decrease r near r c , 
as is clearly seen in Fig. 1. The other important point to notice is that the 
continuum bound ( |l9| ) is satisfied even for a ratio f/b as large as 1/2. This 
indicates that the extra effort required to use the improved action is clearly 
managable in this case. 

2 Conclusions 

We have presented a new action for non-abelian gauge theories which is 
better suited to studying topology on the lattice than the standard Wilson 
action. The idea is to use the Wilson action of an interpolation to a finer 
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lattice of the original lattice configuration. In this way, if the ratio of lattice 
spacings is small enough, the continuum bound on configurations with non- 
zero topological charge is satisfied. In non-abelian gauge theories this is a 
sufficient condition for proper scaling of the topological susceptibility. We 
also considered a new action for the 0(3) cr-model in two dimensions along 
the same lines. A numerical study of discretized continuum instantons in 
this model indicates that already for a ratio of lattice spacings of 1/2, the 
continuum bound is satisfied. Although the results for the 0(3) model are 
very promising, a separate numerical analysis is needed in the Yang-Mills 
case to determine the ratio there. 
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